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Abstract 



A report is presented on Suzaku observations of the ultra-luminous X-ray source X-1 in the starburst 
galaxy M82, made three time in 2005 October for an exposure of ~ 30 ks each. The XIS signals from a 
region of radius 3' around the nucleus defined a 2-10 keV flux of 2.1 x 10^^^ erg s~^ cm~^ attributable to 
C3 \ point sources. The 3.2-10 keV spectrum was slightly more convex than a power-law with a photon index 

of 1.7. In all observations, the HXD also detected signals from M82 up to ^ 20 keV, at a 12-20 keV flux 
of 4.4 X 10"i2 erg s'^ cm'^. The HXD spectrum was steeper than that of the XIS. The XIS and HXD 
^ ■ spectra can be jointly reproduced by a cutoff power-law model, or similar curved models. Of the detected 

' wide-band signals, 1/3 to 2/3 are attributable to X-1, while the remainder to other discrete sources in 

i M82. Regardless of the modeling of these contaminants, the spectrum attributable to X-1 is more curved 

CO ' than a power-law, with a bolometric luminosity of (1.5 — 3) x 10^° erg s^^. These results are interpreted 

, as Comptonized emission from a black hole of 100 — 200 solar masses, radiating roughly at the Eddington 

0^ ' luminosity. 

' Key words: balck hole physics — galaxies:individual:M82 — X-rays:stars 

OO ■ 

o 

1. Introduction where the established mass-accreting BHs, including both 

.,— , , stellar- mass ones and active galactic nuclei (AGNs), arc 

' Ultra- luminous compact X-ray sources (ULXs; e.g., known to emit a considerable fraction of their radiative 

JH '\ Makishima et al. 2000), found in many nearby galax- luminosity. 

ies, are often considered as candidate intermediate-mass Among a number of ULXs known so far, so-called source 

black holes (BHs), because their typical luminosity reach- X-1 in the starburst galaxy M82 is the most extreme ob- 

ing 3 X 10'^^ to 10^^ erg s~^, calculated assuming isotropic ject, as the luminosity reached 10^^ erg s^^. Although the 

radiation, largely exceeds the Eddington limit for known host galaxy M82 has been observed repeatedly with the 

stellar-mass BHs. Since the first discovery by the Einstein past X-ray observatories due to its nearby location (~ 3.6 

satellite (Fabbiano 1989), ULXs have been studied exten- Mpc; Freedman et al. 2004) and its active star forma- 

sively in soft X-rays with ROSAT (Colbert, Mushotzky tion, the bright diffuse thermal X-rays from this galaxy 

1999; Colbert, Ptak 2002; Liu, Bregman 2005; Liu et al. hampered the detection of X-1 till the ASCA era. The 

2006), and then in energies up to ~ 10 keV with ASCA presence of this object, X-1, was first suspected by Tsuru 

(Makishima et al. 2000; Mizuno et al. 2001), Chandra (e.g. et al. (1997), based on differences of fluxes measured with 

Swartz et al. 2004), XMM-Newton (e.g. Feng, Kaaret EXSOAT, Ginga, and ASCA. Subsequently, the object 

2005; Stobbart et al. 2006), and Suzaku (Mizuno et al. was confirmed as a variable hard source through multiple 

2007; Isobe et al. 2008). Searches for their counterparts ASCA observations, and interpreted at first as an active 

have also been carried out in the optical (e.g. PakuU, galactic nucleus veiled by the thermal plasma emission 

Mirioni 2003; Ptak et al. 2006) and radio (e.g. Sanchez- (Matsumoto, Tsuru 1999; Ptak, Griffiths 1999). However, 

Sutil et al. 2006) wavelengths. In contrast, we have lit- thanks to the superb angular resolution of Chandra, it was 

tie information on ULXs in hard X-rays above 10 keV, later revealed to be a non-AGN source located ~ 170 pc 
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off the dynamical center of the galaxy (Matsumoto et al. 
2001; Kaaret et al. 2001; Tsuru et al. 2004). 

To explain the highest luminosity (1 x 10^-'^ erg s~^ in 3- 
20 keV; Rephaeh, Gruber 2002) recorded from M82 X-1, 
a mass of at least 700 Mq is necessary assuming isotropic 
radiation with a sub-Eddington luminosity. Moreover, the 
54 mHz X-ray quasi-periodic oscillation (QPO), detected 
with XMM-Newton and RXTE (Strohmayer, Mushotzky 
2003; Fiorito, Titarchuk 2004; MucciareUi et al. 2006; 
Dewangan et al. 2006; Kaaret et al. 2006), suggests that 
the source, presumably a BH, has a mass of 100-300 Mq, 
assuming the QPO frequency to scale inversely with the 
BH mass. These results, together with a recent study of 
variability pattern (Casella et al. 2008), make M82 X-1 the 
most promising candidate for an intermediate-mass BH. 
It is also suggested to have a 62 day periodicity (Kaaret 
et al. 2006). 

Spectral studies of M82 X-1 have been conducted ex- 
tensively. With ASCA and XMM-Newton, however, such 
attempts were limited to relatively narrow bands, e.g., 
2-10 keV, because of the surrounding bright thermal X- 
rays. As a result, the observed spectra were reproduced 
by a variety of spectral models, including a power-law 
(PL) of photon index F = 1.7 — 2.5 (Matsumoto, Tsuru 
1999; Dewangan et al. 2006), a thermal plasma emission 
model with a temperature 4-11 kcV (Matsumoto, Tsuru 
1999), an unsaturated Comptonization model (Fiorito, 
Titarchuk 2004; Agrawal, Misra 2006), and slim-disk emis- 
sion (Okajima et al. 2006). Since these modelings imply 
vastly different physical conditions of the emission region, 
the nature of M82 X-1 has remained ambiguous. For ex- 
ample, the slim-disk modeling by Okajima et al. (2006) 
leads to a BH mass of ^ 30 Mq, which is significantly 
lower than is implied by the other arguments. Another 
complication with X-1 is that it is heavily confused with 
other point sources in the core region of M82. 

These confusion problems have been resolved by 
Chandra with its high angular resolution. In fact, Kaaret 
et al. (2006) found that the 0.3-8 keV band spectrum of 
X-1 can be represented by a PL model of a photon index 
r 1.7. Nevcrthelss, the information is still insufficient 
to arrive at a clear scenario as to the nature of X-1. 

Evidently, the detection of X-1 in harder energies is 
of crucial importance. The Large Area Counter onboard 
Ginga, sensitive over the 2-30 keV energy band, de- 
tected signals up to ~20 keV from the M82 galaxy, and 
yielded a spectrum which is better described by a thermal 
brcmsstrahlung model of temperature ^5.7 keV than by a 
single PL model (Tsuru 1992). The emission, interpreted 
at that time as arising from a thin hot plasma, is now 
more likely to be attributable primarily to X-1. Cappi et 
al. (1999) successfully described a 3-30 keV BeppoSAX 
spectrum of M82 by a thermal plasma model of tem- 
perature ~ 8.2 keV. Furthermore, RXTE observations of 
M82 on 30 occasions yielded 3-50 kcV spectra (Rephaeli, 
Gruber 2002); the spectra acquired in brighter and fainter 
phases were both reproduced by thermal plasma models, 
with a temperature of ~ 6.6 and 7.4 keV, respectively. 
These existing results in harder energies commonly in- 



dicate mildly curving broad-band spectra, rather than 
a PL-like one. However, the fields-of-view (FOVs) of 
these non-imaging instruments also contained the spiral 
galaxy M81, hosting a low-luminosity AGN (Ishisaki et al. 
1996; Pellegrini et al. 2000; La Parola et al. 2004), located 
37' away from M82. Therefore, possible contamination to 
the observed signal was not excluded. 

The Hard X-ray Detector (HXD; Takahashi et al. 
2007; Kokubun et al. 2007) onboard Suzaku (Mitsuda et 
al. 2007), with its unprecedented sensitivity in a broad 
(10-600 keV) energy band and a tightly coUimated FOV 
of 34' X 34' (FWHM), is expected to provide improved 
hard X-ray information on M82 X-1. Here, we report on 
the Suzaku detection of M82 X-1 up to ^ 20 keV, and in- 
terpret the results in terms of an intermediate-mass BH. 

2. Suzaku Observations and Data Reduction 

We observed M82 with Suzaku three times in the 
Science Working Group (SWG) phase; on 2005 October 
4 UT11:50 through October 5 UT03:24, October 19 
UT00:40~22:32, and October 27 UT11:05 through 28 
UT02:24. Aiming primarily at the study of hot-gas out- 
flows from M82, the three observations were all con- 
ducted with the XIS optical axis pointed at {a^°°° ,5"^°°°) = 
(9'^55™34%69°45'53"), which is ~ 5' north of the galaxy 
The XIS results on this offset region have already been 
pubhshed as Tsuru et al. (2007). In the present paper, 
we focus on the central region of this galaxy, including in 
particular X-1. 

Figure 1 shows an XIS image in the 2-10 keV band 
from the second observation. The X-ray bright part coin- 
cides with the central region of M82, including luminous 
X-ray sources such as X-1. The M81 nucleus during this 
observation was right on the edge of the HXD-PIN FOV, 
where the transmission is < 1% even considering attitude 
uncertainties. Similarly, the HXD-PIN transmission for 
M81 was < 3% both in the first and third observations. 
As shown in the figure, two ULXs known in the present 
sky region, M81 X-6 and Holmberg XI X-1 (= M81 X-9; 
Fabbiano 1988a), were both outside the HXD-PIN FOV. 

We used the software package HEADAS version 6.3.1 
and xselect version 2.4 to extract light curves, images, 
and spectra, and xspec version 11 to analyze the spectra. 
The XIS and HXD datasets were screened by the revi- 
sion 1.2 pipeline processing using the standard criteria as 
follows. Events were discarded if they were acquired un- 
der the South Atlantic Anomaly (SAA), or at regions of 
low cutoff rigidity (< 6 GV for XIS and < 8 GeV c'^ for 
HXD), or with the Earth elevation angle of < 5°, or dur- 
ing those periods when the telemetry was saturated. The 
XIS events were further excluded when the source eleva- 
tion above the sunlit Earth is < 20°, and finally events 
with grade 0, 2, 3, 4, and 6 were selected. After these 
screenings, the three observations gave effective exposure 
of 29 ks, 36 ks, and 26 ks with the XIS, while 24 ks, 31 
ks, and 22 ks with the HXD. 
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3. Data Analysis and Results 

3.1. XIS Data Analysis 

We extracted on-source XIS events from a region of 
3' radius centered on the M82 nucleus, and background 
events from another region of 3' radius which is free from 
the diffuse X-ray emission associated with M82. In figure 
1, these two regions are indicated by dashed circles. We 
utilize events from all four CCD cameras, after co-adding 
those from the three FI cameras (XIS 0, XIS 2, and XIS 3) 
which have virtually identical responses. 

3.1.1. XIS light curves 

Figure 2 shows 2-10 kcV XIS light curves obtained 
in the three observations. The background counts are 
negligible (< 2%), compared to the signal counts which 
were determined as 0.460 ± 0.007, 0.489 ± 0.006, and 
0.425 ±0.007 cts s~^ in the three observations (the er- 
rors being 90% confidence ranges). Thus, the signal on 
the second and third occasions arc higher by ~ 6% and 
lower by ~ 8%, respectively, than that in the first observa- 
tion. By fitting these light curves with a constant, we ob- 
tained xVd-O.f (degree of freedom) of 95.7/95, 190.0/131, 
and 110.1/87, respectively. The probability of these signal 
counts being consistent with a constant value is ~ 40%, 
1%, and 5%, in the first, second, and third observations, 
respectively. Therefore, the source is likely to have var- 
ied during the second observation; in fact, the light curve 
suggests a gradual intensity increase by ~ 10% across the 
gross coverage of about one day. 

3.1.2. XIS spectra 

Figure 3 shows 0.5-10 keV spectra from the FI CCD 
(black; averaged over XIS , XIS 2, and XIS 3) and the BI 
CCD (gray; XIS 1), obtained by accumulating events over 
the source region of figure 1, and then subtracting those 
from the background region. To grasp the rough spectral 
properties, we tentatively summed up the three observa- 
tions. Thus, the spectra are dominated by a rather hard 
continuum at energies above ~ 3 keV, whereas by line-rich 
softer thermal emission at lower energies. As confirmed 
through past observations, the hard continuum must be 
coming from luminous point X-ray sources including X- 
1 in particular, while the soft thermal emission from the 
diffuse hot plasmas (e.g., Tsuru et al. 1997; Moran et al. 
1997; Griffiths et al. 2000; Stevens 2003; RanaUi et al. 
2008; Strickland, Heckman 2007). A rough evaluation us- 
ing a thermal model and a PL continuum indicates that 
the XIS spectra in energies above 3.2 keV (i.e., above Ka 
line at 3.14 keV from H-likc Ar) are contaminated by the 
thermal component by no more than 9%. Since we are 
interested in X-1, below we study the XIS spectra in the 
3.2-10 keV range, discarding their softer part. 

In the hard-band XIS spectra, we observe an emission 
line structure at 6-7 keV; this must be Fe-K lines, as re- 
ported in previous observations (e.g. Tsuru et al. 1997; 
Matsumoto, Tsuru 1999; Cappi et al. 1999; Rephaeli, 
Gruber 2002; Strohmaycr, Mushotzky 2003; Agrawal, 
Misra 2006). Using the Chandra and XMM-Newton data, 
Strickland, Heckman (2007) spectrally resolved the iron 
line complex into 6.7 keV (from He-like iron) and 6.4 keV 



(from nearly neutral iron) components, and showed that 
both are primarily distributed diffusely rather than asso- 
ciated with discrete sources. 

3.2. HXD- PIN Data Analysis 

Through the method described in section 2, we ex- 
tracted the HXD-PIN data from the three observations. 
Although the data covers an energy range of 10-70 keV, 
we hereafter discard the lowest energy band of 10-12 keV, 
where the data are often affected by increased noise. Light 
curves and spectra are corrected for instrumental dead 
times after Takahashi et al. (2007) and Kokubun et al. 
(2007). In the present paper, we do not intend to analyze 
the HXD-GSO data, since the signals from M82 are esti- 
mated to be below the GSO detection limit. In addition, 
the HXD-GSO FOV becomes progressively wider above 
~ 80 kcV, and hence the contamination by M81 becomes 
difficult to eliminate. 
3.2.1. Background subtraction 

To perform background subtraction which is crucial 
for the HXD analysis, we utilized non X-ray background 
(NXB) models provided by the HXD detector team 
(Kokubun et al. 2007). When the present data anal- 
ysis was carried out (Miyawaki 2008), there were two 
NXB models available, so-called bgd_A (Watanabe et al., 
suzakumemo-2007-01) and bgd_D (dead-time uncorrected; 
Fukazawa et al., suzakumemo-2007-02), which employ dif- 
ferent algorithms. We produced fake background events 
using the two NXB models, both under the same observ- 
ing conditions as the actual data, and both utilizing the 
same auxiliary information such as the measured upper- 
discriminator hit rates. By processing the fake events 
through the same screening criteria as for the actual on- 
source data, we derived two NXB event sets for each ob- 
servation. Then, bgd-A and bgd-D respectively gave 12-70 
keV count rates of 0.46 and 0.45 cts s~^ on October 4, 0.49 
and 0.48 cts s"^ on October 19, and 0.52 and 0.50 cts s"^ 
on October 27. Thus, the two NXB models agree within 
4%. Hereafter, we employ bgd_A for our spectral anal- 
ysis, because it tends to systematically over-predict the 
NXB than the other model in this particular observation, 
and hence is expected to yield more conservative results 
for sources hke M82 X-1 that are close to the HXD-PIN 
detection limit. 

We examined the NXB modeling accuracy utilizing 
Earth-occulted data, assuming the Earth to be "dark" 
for HXD-PIN. These data were extracted using the same 
criteria as for the on-source data, except that the Earth 
elevation angle was required to be < —5° instead of > 5°. 
The attained exposure was 5.3 ks on October 19 and 1.3 
ks on October 27, while the October 4 data were free 
from Earth occultations. Figure 4 compares the observed 
Earth-occultation data, with the synthetic bgd-A spectra 
which emulate the NXB to be observed during these oc- 
cultation periods. The observed and modeled NXB count 
rates in the 12-70 keV energy range are 0.56 ±0.02 and 
0.55 cts s~^ respectively on October 19, and 0.54 ±0.04 
and 0.54 cts s^^ on October 27 likewise, where the errors 
refer to 90% confidence levels. Statistical errors associated 
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with the NXB model are neglected, because the model 
produces 10 times larger number of fake events. When 
the energy band is limited to 12-20 keV, where signals 
are detected as described later, the observed and modeled 
NXB count rates become 0.31 ±0.01 and 0.31 cts s^^ re- 
spectively on October 19, and 0.30 ±0.03 and 0.30 cts s~^ 
on October 27 likewise. 

As confirmed above, the Earth-occultation data and 
the modeled NXB data agree well with each other within 
statistic errors, which are 4% and 10%, on October 19 and 
27, respectively. This statement applies both to the 12-70 
keV and 12-20 keV bands. Further considering the NXB 
modeling accuracy of HXD-PIN which is generally ^ ±5% 
(Mizuno et al. 2007), and the 4% difference between the 
two NXB models, we hereafter regard that the NXB model 
correctly reproduces the HXD-PIN NXB within a system- 
atic error of at most 4%. 
3.2.2. HXD signal count rates 

After the above preparations, we accumulated the on- 
source HXD-PIN data of M82, and derived 12-20 keV 
count rates from the three observations. In table 1, these 
results are compared with the NXB model predictions in 
the same energy band. In all three observations, the on- 
source data exhibit an excess count rate by ~ 0.04 cts s^^ 
above the modeled NXB. Of these excess counts, we ex- 
pect 0.018 cts s~^ to be coming from the cosmic X-ray 
background (CXB), as calculated from the HEAO-1 re- 
sults (Boldt 1987; Gruber et al. 1999). Further subtract- 
ing this CXB contribution, we are finally left with 0.016- 
0.027 cts s"^ (table 1) in 12-20 keV to be attributable to 
the net signals from M82. Ahhough the 12-20 keV CXB 
intensity within the HXD-PIN FOV is expected to fluc- 
tuate by ~ ±12% (rms) from sky to sky, this is only to 
±1% of the NXB in the same energy range, because the 
average CXB count rate therein, typically 7% of the NXB, 
is already low. Thus, the CXB fluctuation is negligible. 

The PIN signals, remaining after subtracting the NXB 
and CXB, have statistical significance of 3 to 5 sigmas 
on individual observations, and 8 sigmas when the three 
data sets are summed together. Furthermore, the implied 
12-20 keV signal rates, ~ 0.02 cts s-\ amount to 6-10% 
of the NXB count rate in the same band. This exceeds 
the typical NXB reproducibility during the present obser- 
vations, 4%, as estimated in section 3.2.1. Therefore, the 
positive signals are thought to be significant even consid- 
ering the systematic NXB errors. 

After we had completed our data analysis, the back- 
ground model bgd_D was updated, by taking into ac- 
count dead-time correction and improving the esti- 
mation method (Mizuno et al., suzakumemo-2008-03). 
Accordingly, we re-generated the background files using 
this improved method, to find that it gives essentially the 
same results as the previous hgd_D. Therefore, we retain 
our original choice of bgd-A, which tends to overestimate 
the NXB and hence gives more conservative results. We 
later examine the effects of these different background 
models by artificially changing the NXB by ±4%. 



3.2.3. HXD light curves 

As described in section 3.1, the XIS signal from the 
central region of M82 exhibited rather small (if any) vari- 
ations during the present observations. Supposing that 
the positive HXD-PIN signals, as revealed above, are also 
coming from an approximately identical region in M82, 
we expect rather constant count rates as well. If, on the 
other hand, the positive HXD-PIN counts were an artifact 
of wrong NXB subtraction, they would vary with time 
typically on time scales of minutes to hours, just refiect- 
ing those in the NXB (Kokubun et al. 2007). For this 
purpose, the 12-20 keV PIN light curves of the three ob- 
servation were derived, and presented in figure 5 together 
with the behavior of the modeled NXB. 

On all occasions, the background-subtracted counts in 
figure 5 thus show on-average positive values, which agree 
with those in table 1. Furthermore, the net signal (green) 
appears relatively constant within statistical errors, in 
contrast to significant raw-count variations (black) that 
reflect those in the NXB (red). When fitted with a 
constant, these background-subtracted light curves yield 
xVd.o.f of 41.7/28, 42.8/39, and 31.2/30, in the first, 
second, and third observations, respectively. These imply 
probabilities of ~ 5%, 35%, and 40% for the fight curves to 
be consistent with being constant. Thus, we consider that 
the positive HXD-PIN counts listed in table 1 persisted 
throughout the observations, while the variable NXB was 
removed successfully. 

In this way, we have carefully examined the HXD-PIN 
signals using the Earth-occultation data and the on-source 
light curves. The results confirm that the positive HXD- 
PIN signals above the CXB are real, and are coming from 
the M82 sky region. 

3.2.4. HXD spectra 

Figure 6 shows the full-range PIN spectra (black) from 
the three observations, compared with the corresponding 
NXB models {bgd_A; red). The figure also shows the NXB- 
subtracted data (green), and the predicted CXB spectrum 
(blue) based on the HEAO-1 results (Boldt 1987; Gruber 
et al. 1999). Thus, after the NXB subtraction, positive 
signals remain up to ^ 30 keV in all three observations. 
Furthermore, the signal exceeds the expected CXB at least 
up to ~ 20 keV. In the 12-20 keV band, these spectra 
imply an average signal count rate of ^ 0.023 cts s"-'^ above 
the CXB. This value agrees with that derived from the 
simple count-rate measurements (section 3.2.2; table 1), 
and the light curve analysis (section 3.2.3). 

3.3. Spectral Model Fitting Analyses 

In this subsection, we perform model fitting first to the 
XIS spectra, and then to those of HXD-PIN, both aver- 
aged over the 3 days. Afterwards, the XIS and HXD-PIN 
spectra are jointly analyzed. 
3.3.1. Model fits to the XIS spectra 

Among the three observations, the 3.2-10 keV XIS spec- 
tra differ to some extent in their normalization; as inferred 
from figure 2, those of the 2nd and 3rd observations are 
higher by 2.1% and lower by 10.6%, respectively, than 
that of the first occasion. However, after scaling by these 
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factors, the three spectra agree with one another within 
their errors, which are typically 5% per channel under the 
present binning. Therefore, we hereafter analyze the XIS 
spectra summed over the 3 observations. 

We jointly fitted the 3.2-10 keV XIS FI and BI spectra, 
thus summed over the three observations, with a model 
consisting of a PL for the continuum and two Gaussians 
for the Fc-Kq, lines. The XIS response file and the asso- 
ciated ancillary file was generated with "xisrmfgen" and 
"xissimarfgen" tool, respectively (Ishisaki et al. 2007). 
We constrained the overall model normalization to be the 
same between the FI and BI detectors. The hydrogen col- 
umn density describing photoelectric absorption is here 
and hereafter fixed at 1.1 x 10^^ cm~^, as determined with 
Chandra (Kaaret et al. 2006) under the least contamina- 
tion by the diffuse emission: this amount of absorption 
attenuates the continuum at 3.2 keV only by 10%, and 
if treated as a free parameter, the absorption would cou- 
ple too strongly with the continuum slope. The model 
yielded a fit as shown in figure 7 (left), with the best-fit 
parameters summarized in table 2 under the column of 
"XlS-only". Thus, the PL photon index became P = 1.67. 
Although not given in the table, one Gaussian is cen- 
tered at an energy of 6.67 ±0.02 keV, with an intrinsic 
width (in Gaussian sigma) of < 0.11 keV and an equiv- 
alent width (EW) of 59± 11 eV. Its inclusion improved 
the fit by Ax^ ~ 88. The other Gaussian, with its energy 
fixed at 6.4 keV and assumed to be narrow, exhibited an 
EW of 19 ± 8 eV, and improved the fit by Ax^ ~ 15. 

The 2-10 keV flux detected with the XIS, averaged over 
the three data sets, becomes 2.1 x 10~^^ erg cm~^, 
which yields a 2-10 keV luminosity of 3.4 x 10*° erg s~^ 
at 3.6 Mpc. This is typical of, or relatively lower than, 
those measured previously from narrow regions (depend- 
ing on the angular resolution) around X-1, which scat- 
ter over (1 — 6) X 10~^^ erg s~^ cm~^ (e.g. Griffiths et 
al. 1979; Watson et al. 1984; Fabbiano 1988b; Shcaaf 
et al. 1989; Tsuru 1992; Tsuru et al. 1997; Matsumoto, 
Tsuru 1999; Ptak, Grifiiths 1999; Cappi et al. 1999; 
Matsumoto et al. 2001; Kaaret et al. 2001; Rephaeli, 
Gruber 2002; Strohmayer, Mushotzky 2003; Kaaret et al. 
2006; Dewangan et al. 2006). 

Although the above model approximately reproduces 
the observed XIS spectra, the fit is formally unacceptable 
with x^/i^ = 1.83 (table 2). As shown in figure 7a, the 
model overpredicts the data above ~ 7 keV, suggesting 
the spectrum to have a more convex shape than a sim- 
ple PL. Actually, the best-fit photon index becomes dis- 
crepant as 1.56 ±0.03 and 2.0 ±0.2, when the fit energy 
range is restricted to 3.2-6.0 keV and 7-10 keV, respec- 
tively. The PL fit could be made acceptable by allowing 
the low-energy absorption to vary freely, but then the col- 
umn density must become ^ 3.7 x 10^^ cm~^, much ex- 
ceeding the Chandra value of 1.1 x 10^^ cm^^ (Kaaret et 
al. 2006). Therefore, the intrinsic continuum is suggested 
to be curved. To express such a convex shape, we fixed 
again the absorption to 1.1 x 10^^ cm~^, and replaced the 
PL continuum by a cutoff PL model, namely a PL multi- 
plied by an exponential cutoff factor. Then, as shown in 



figure 7b and described in table 2, the fit has become ac- 
ceptable with x^/'^ ~ 1-06, yielding P ~ 0.76 and a cutoff 
energy of Ecut = 5.8 keV. Therefore, we conclude that the 
spectrum has a more convex shape than a PL. 

3.3.2. Model fits to the HXD-PIN spectra 

We conducted similar model-fitting analyses to the 
background (NXB±CXB) subtracted 12-20 keV PIN 
spectra, again summed over the three observations. The 
fit utilized a response matrix at the HXD nominal po- 
sition, ae_hxd_pinhxnom_20060814.rsp, and an ancillary 
response file generated by the hxdarfgcn tool. The latter 
takes into account the fact that the HXD FOV center is 
~ 5'. 7 offset from the M82 nucleus. Then, a PL model 
gave a successful fit with x^/'^ = 9.2/6, yielding a pho- 
ton index of P = 3.4t|!) g (statistical 90% confidence er- 
rors), and a 12-20 keV fiux of 4.4 x 10"^^ erg s"^ cm"^. 
These become P = S.lttl and 1.9 x lO^^^ erg s^^ cm"^ 
if the NXB is artificially increased (over-subtracted) by 
4%, while P = 3.ltoi and 6.7 x lO'^^ ^^.g g-i ^^-2 ^^^^ 
it is reduced (under-subtracted) by 4%. In any case, the 
implied spectral slope (P = 3 — 4) is significantly steeper 
than that found in the XIS range (P = 1.67), reinforcing 
the gradual spectral steepening found in the XIS band. 

3.3.3. Joint fits to the XIS and PIN spectra 

In order to best utilize the wide-band capability of 
Suzaku, we fitted the XIS and HXD spectra simultane- 
ously, using respective averages over the three observa- 
tions. Since the energy range becomes much expanded, we 
tried a larger variety of continuum models than employed 
so far. In all modelings, the two Gaussians were included 
in the same manner as in section 3.3.1. The model normal- 
ization was constrained to be the same between the XIS 
and HXD-PIN, because the HXD-PIN ancillary response 
file (section 3.3.2) has already taken into account the dif- 
ferent FOVs between the two instruments (assuming that 
the HXD-PIN signal is coming from the M82 central re- 
gion). The obtained fits are shown in figure 8, while the 
fit parameters are summarized in table 2 under the col- 
umn of "XIS±PIN" . (The table also gives the results ob- 
tained without using the HXD data.) Furthermore, like 
in section 3.3.2, we considered the case where the NXB of 
HXD-PIN is systematically changed by ±4%. 

First, we tried a PL model, but the joint fit was not 
acceptable with x^/^^ ^ 2.6 even allowing for the 4% 
NXB uncertainty (table 2). As shown in figure 8a, the 
model over-predicts the XIS and HXD-PIN data in the 
energy ranges of > 7 keV. Although the fit is improved to 
X^/i' 1.9 by tentatively allowing the model normaliza- 
tion to take separate values between the two instruments, 
the HXD vs. XIS normalization ratio then became < 0.7. 
This is too small, considering the XIS vs. HXD calibra- 
tion uncertainties which are at most ~ 15% (Kokubun et 
al. 2007). Thus, the joint fit reinforces the inference made 
so far, that the X-1 spectrum is more convex than a PL. 
Below, we return to our original procedure of equalizing 
the model normalization between the two instruments. 

Next, we tried an empirical thermal bremsstrahlung 
(bremss) model. Although the fit is better jv = 
1.3 — 2.2 depending on the NXB systematic errors) than 
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the PL fit, the model is still over-predicting in the > 7 keV 
energy range as shown in figure 8b. The obtained temper- 
ature, T ~ 14 keV, is higher than the values of 5 — 8 keV 
measured with Ginga (Tsuru 1992), BeppoSAX (Cappi 
et al. 1999), and RXTE (RephaeH, Gruber 2002). These 
results remain essentially unchanged even when changing 
the PIN NXB by ±4%. 

Third, we employed a cutoff PL model. As shown in 
figure 8c and table 2, it has given jv = 1.07— 1.26, to- 
gether with F = 0.5 - 0.9 and E^^t = 4-8 keV. Therefore, 
this model remains successful even including the HXD- 
PIN data. Because a cutoff PL gives an approxima- 
tion to an unsaturated thermal Comptonization process 
of some soft seed photons, we replaced it with a ther- 
mal Comptonization model, called "CompTT" (Titarchuk 
1994) in xspec. Since the input soft-photon energy is 
not well constrained by the present data, we fixed it to 
several values between 0.1 and 0.3 keV. Then, as pre- 
sented in figure 8d, the fit was similarly acceptable with 

jv = 1.09— 1.32, and gave an electron temperature of 
To = 2.3 — 2.9 keV and an optical depth of t = 7 — 9; the re- 
sults were essentially independent of the soft-photon tem- 
perature in the above-mentioned range. The obtained val- 
ues of Te and r are both consistent with the XMM-Newton 
results on M82 X-1 by Agrawal, Misra (2006). 

Finally, to evaluate a scenario of emission from an 
optically-thick accretion disk, we tested a "disk black- 
body" (diskbb) model which represents multi-color black- 
body emission from a standard accretion disk (Mitsuda et 
al. 1984; Makishima et al. 1986). As shown in figure 8e, 
the fit was rather poor with x^/i^= 1.32 — 1.72, because the 
model has too convex a shape and hence under-predicts 
the data in the higher energy range. To improve the 
fit, we then generalized the diskbb model to "variable-p" 
model (Mineshige et al. 1994; Hirano et al. 1995; Kubota, 
Makishima 2004; Kubota et al. 2005; Okajima et al. 
2006; Mizuno et al. 2007). It assumes that the disk tem- 
perature scales as , where r is the radius and p is a 
free parameter, with p = 0.75 corresponding to diskbb. It 
gives a good approximation to emission from a slim disk 
(Watarai et al. 2000; Ebisawa et al. 2003), which is ex- 
pected to form under very high accretion rates. The fit 
has then become acceptable with /v=l .09 — 1 .24 (figure 
8f), and the obtained parameters, namely the inner-disk 
temperature Tin = 3.4 — 3.6 keV and p = 0.61 — 0.65, are 
consistent with those of Okajima et al. (2006) who ap- 
plied the the same model to the XMM-Newton data of 
X-1 acquired at a similar luminosity to the present case. 

Ftom these studies, we have found three models that 
successfully and simultaneously reproduce the XIS and 
HXD-PIN data. A cutoff PL model with E^nt ~ 5 keV, a 
CompTT model invoking a relatively cool and thick elec- 
tron cloud, and a variable-p model with a very high disk 
temperature. In all these cases, the two Gaussians were 
found to take nearly the same parameters as found with 
the XIS (section 3.3.1). 

As seen from table 2, the addition of the HXD-PIN 
data have little affected the best-fit parameters of these 
successful models, with similar or somewhat reduced er- 



rors. In contrast, their inclusion has allowed us to much 
more securely exclude the other unsuccessful models. In 
this sense, the XIS and HXD-PIN data are fully consis- 
tent with each other. When the NXB level is artificially 
changed by ±4%, the fit goodness of all the three suc- 
cessful models degrades considerably. This reinforces the 
correctness of the NXB subtraction, unless the detected 
signal has a very complex broad-band continuum shape 
which requires a composite modeling. 

3-4-. Estimation of Contamination from Other Sources 

Although we have so far quantified the broad-band 
Suzaku spectra, the results may not be totally ascribed 
to X-1, because the HXD-PIN data are likely to be con- 
taminated by M81, as well as by other discrete sources 
in M82. The latter must contaminate the XIS data as 
well, because many of them are clustered in too narrow 
a region, ~ 30" in radius around the M82 center, to be 
resolved with the XIS. 

3.4.1. M81 

As described in section 2 and illustrated in firgure 1, 
the M81 nucleus was on the edge of the HXD-PIN FOV, 
with transmission of at most 1 — 3%. Therefore, the con- 
tamination is expected to be small, but may not be totally 
neglected. 

Using BeppoSAX, Pellegrini et al. (2000) measured the 
10-100 keV flux of the M81 nucleus as 7.4 x 10"" erg s"! 
cm~^, and described its spectrum by a PL with F ^ 1.84. 
The 12-20 keV flux is impfied to be 1.5 x 10"" erg s^^ 
cm~^. According to La Parola et al. (2004), furthermore, 
the M81 nucleus varied in the 0.5-2.4 keV band by a factor 
of 4 through past observations, and became up to ^ 1.4 
times as bright as observed by Pellegrini et al. (2000). 

Adopting an extreme assumption that the M81 nucleus 
was at its brightest level during the present observations, 
and that the spectrum retained F ~ 1.8, the angular trans- 
mission described in section 2 predicts this low-luminosity 
AGN to contribute up to < (2 - 6) x 10"" erg s'^ cm'^ 
to the 12-20 keV HXD-PIN data. These correspond 
to 5 — 15% of the background-subtracted PIN flux (sec- 
tion 3.3.2). Taking an exposure-weighted mean over the 
three occasions, the M81 contamination to the HXD-PIN 
signal is estimated to be at most 10%. Since the PIN 
spectra in figure 8 have statistical errors of several tens 
percent per bin, and since M81 may not have been in the 
brightest phase, we conclude that the effect of the M81 
nucleus is negligible. This is an advantage of the HXD 
with its tightly coUimatcd FOV. 

3.4.2. Other discrete X-ray sources in M82 

Besides X-1, the central region within ~ 30" of the 
M82 nucleus harbors about a dozen X-ray point sources 
more luminous than ~ 10^^ erg s~^ (Griffiths et al. 2000; 
Matsumoto et al. 2001; Kaaret et al. 2006; Strickland, 
Heckman 2007). To study their fluxes and spectra, and es- 
timate their contamination to the present Suzaku results, 
we re- visited the Chandra archival data. By 2007 August, 
the central regions of M82 were observed six times with 
the Chandra ACIS for exposures longer than 10 ks; these 
are ObsID 361, 1302, 2933, 6097, 5644, and 6361. Using 
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the CIAO 3.4 software package and standard screening cri- 
teria, we extracted point sources which are located within 
2'. 9 of X-1; these can contaminate both the present XIS 
and HXD spectra. Referring to the ROSAT PSPC image, 
any X-ray sources that are outside the XIS integration 
radius but within the HXD-PIN FOV arc considered neg- 
ligible. 

Through this analysis, we confirmed the report by 
Kaaret et al. (2006) that a second brightest source, which 
they call X42. 3-1-59, is located at ~ 5" off X-1; we here- 
after call it X-2. Across the six Chandra observations, 
the 2-10 kcV flux of X-2 varied largely, from less than 
~ 10~^^ erg s~-^ cm~^ to 6.9 x 10"^^ erg s~^ cm~^, with 
the latter amounting to 1/3 of the 2-10 keV XIS signal 
flux. As already reported by Kaaret et al. (2006), the 
3-8 keV spectrum of X-2 at its brightest phase was de- 
scribed by a PL with F 1.3, absorbed by a column of 
A^H = (2 — 3) X 10^^ cm~^, and a soft-excess component 
was required if lower energies are included. 

Fainter X-ray sources, other than X-1 and X-2, were 
found to have a 2-10 keV luminositiy of ^ 3 x 10^^ erg 
s~^ each. Among the six Chandra observations, their 
summed 2-10 keV flux remained relatively constant at 
(7 — 9) X 10~^^ erg s^^ cm~^, which accounts for about 
1/3 of the XIS signal. Their summed 3-8 keV spectrum 
was reproduced by a PL with F = 1.5 — 1.7. 

As shown so far, we expect the fainter discrete sources 
to account for about 1/3 of the detected XIS signal. 
Furthermore, another 1/3 could be ascribed to X-2 if it is 
in the brightest phase. Therefore, the soft X-ray flux ac- 
tually attributable to X-1 is estimated to be 2/3 (if X-2 is 
very faint) to 1/3 (if X-2 is brightest) of the detected 2-10 
keVXIS signal of 2.1x10"" erg s"^ cm'^ (section 3.3.1). 
These contaminants are expected to introduce a larger 
uncertainty to the HXD-PIN results, because their hard 
X-ray spectra cannot be measured separately. Below, we 
consider this issue assuming two extreme cases. 
3.4-3. Re-evaluation of the spectra considering contami- 
nation 

We first consider an extreme case where X-2 is faint 
enough to be neglected, and hence we need to consider 
only the fainter sources. If their summed spectrum de- 
termined with Chandra (a PL with F = 1.5 — 1.7) were 
extending into the hard X-ray range, the implied counts 
at 30 keV would reach - 10% of the HXD-PIN back- 
ground, which is well above the systematic background 
uncertainty. Furthermore, the 30- 40 keV signals, at- 
tainable with HXD-PIN in the 77 ks exposure (the three 
observations summed), would have a statistical signifi- 
cance exceeding 6 sigmas. Thus, we must have detected 
their signals at least up to ~ 40 keV, as illustrated in fig- 
ure 9a. However, the actual 12-20 keV PIN spectrum was 
rather soft with F ~ 2.5, and was undetectable above ^ 20 
keV. Therefore, the summed spectrum from the fainter 
sources should bend in energies above ^ 7 kcV. This is 
reasonable, because these sources are considered to be 
mainly luminous low-mass X-ray binaries, of which the 2- 
20 keV spectrum is empirically approximated by a thermal 
bremsstrahlung model with a temperature of T = 7 — 13 



keV (Makishima et al. 1989). Such a bremsstrahlung 
model has a mildly curved shape, and its 3-8 keV portion 
can be approximated by a PL with F 1.5 (Makishima ct 
al. 1989), in agreement with the Chandra measurement. 

With the above consideration, we modeled the summed 
contribution of the faint sources by a T = 10 keV ther- 
mal bremsstrahlung model, of which the 2-10 keV flux 
was fixed at 8 x 10"^'^ erg s~^ cm~^ as specified by the 
Chandra data. Including this model as a fixed com- 
ponent, we jointly re-fitted the same Suzaku XIS and 
PIN spectra using the same family of models as in sec- 
tion 3.3.3. As shown in tabic 3 under the column "Case 
1" , and presented in flgure 9b, the obtained results are 
close to those obtained before. However, the inferred 
X-1 spectrum became softer, because the T = 10 keV 
thermal bremsstrahlung model, which is now subtracted 
away, is somewhat harder than the original wide-band 
Suzaku spectrum. As a result, not only the successful 
three models, but also the simple multi-color disk model 
with Tin ~ 3.1 keV, have become acceptable. Equivalently, 
the variable-p model gave p ~ 0.75, which just implies a 
standard accretion disk. After subtracting the contam- 
ination, the 2-10 kcV flux attributable to X-1 becomes 
1.2 X 10"^^ erg s^^ cm^^, and the luminosity therein is 
1.9 X 10"° erg a'^. 

As the other extreme, we may consider the case where 
the present observations happened to catch X-2 at its 
brightest level ever recorded by Chandra, and hence its 
contribution to the Suzaku spectra must be considered as 
well. Although the X-2 spectrum obtained by Chandra 
at the brightest epoch was rather hard (F ~ 1.3), it can- 
not keep extending into the > 10 keV energy range, for 
the same reason as discussed just above for the fainter 
sources. Thus, the X-2 spectrum is inferred to have a 
power-law shape in the 3-8 keV range, together with a 
high-energy turn over, and a soft excess (Kaaret et al. 
2006). Such a spectral composition is typical of "PL 
state" ULXs (e.g., Feng, Kaaret 2005; Stobbart et al. 
2006). Although these objects generally have a steeper 
range of spectra with F = 1.6 — 2.5, a case with F ^ 1.4 
is known (Kubota et al. 2001b). Regarding X-2 as such 
an object, we decided to express its spectrum by a cutoff 
PL model with F = 0.5 and a cutoff energy of i?cut = 5.0 
keV, absorbed with iVn = 2.5 x 10^^ cm~^ (Kaaret ct al. 
2006). The adopted cutoff energy is a typical value found 
among power-law type ULXs (Miyawaki 2008; Mizuno et 
al. 2007), and the photon index of 0.5 was chosen so that 
the model, of the form oc cxp(— i?/£'cut)£'~'", mimics the 
F = 1.3 PL in the 3-8 keV band. Its normalization was 
adjusted to reproduce the Chandra flux. 

Including these two flxed models (one for X-2 and the 
other for the fainter sources), we repeated the fltting pro- 
cedure; the obtained results are presented in flgure 9c, and 
the parameters are given in table 3 under the column de- 
noted "Case 2". Although the 2-10 keV flux attributable 
to X-1 further decreased to 7.0 x 10~^^ erg s~^ cm~^, and 
its 2-10 kcV luminosity to 1.1 x 10^" erg s~^, the spectral 
results have remained mostly the same as the case of X-2 
being neghgible. This is because the assumed X-2 spec- 
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trum is relatively similar in shape to the original Suzaku 
XIS plus HXD-PIN spectra, as is clear from the cutoff PL 
fit parameters in table 2. 

These studies confirm that the removal of the estimated 
contributions from X-2 (under the most extreme assump- 
tion) and the other fainter sources leaves significant 3-20 
keV counts that are attributable to X-1. These remaining 
hard X-ray photons still form a clearly curved spectrum, 
which excludes a PL and a thermal bremsstrahlung mod- 
eling, similar to the case of ignoring the contamination. 
Thus, except for the normalization, our basic results on X- 
1 are relatively unaffected by the contaminating sources, 
in spite of apparently large uncertainties that they intro- 
duce. This is because we would encounter an unphysical 
condition of negative X-1 flux at ~ 20 keV if the con- 
taminants had significantly harder spectra, while we can 
simply utilize the parameters in table 2 if the spectra of 
the contaminants were much softer and hence negligible 
in the HXD-PIN range. 

4. Discussion 

Summary of the results 

With the Suzaku XIS, we detected a 2-10 keV flux of 
2.1 X 10""'^^ erg s~^ cm~^ that is attributable to the as- 
sembly of luminous point sources within 3' of the nucleus 
of M82. Of this flux, about 2/3 can be accounted for 
by the sum of the two most luminous sources, X-1 and 
X-2, whereas the remaining 1/3 is ascribed to the other 
fainter ones. The share of 2/3 is totally attributable to 
X-1 if the highly variable X-2 was negligibly dim at that 
time, whereas it is further subdivided roughly equally be- 
tween X-1 and X-2 if X-2 was at the brightest state ever 
recorded. We hence estimate the 2-10 keV flux of X-1 as 
(7— 14) X 10^^^ erg s^^ cm^^, and its 2-10 luminosity as 
(1.1 - 2.1) X 10"*" erg s^^ at 3.6 Mpc, both with a typical 
systematic uncertainty by a factor of two. The estimated 
luminosity of X-1 is similar to, or somewhat lower than, 
those so far measured from this source. 

From a wider sky region including M82, we also de- 
tected significant hard X-ray emission with the Suzaku 
HXD-PIN, at a 12-20 keV fluxof 4.4x lO'^^ ergs"! cm^^. 
When summed over the three observations, the HXD-PIN 
signal (after subtracting the NXB and CXB) has a statis- 
tical significance of 8 sigmas, and about twice as high as 
the systematic background uncertainty. Therefore, the 
signal is considered real. The HXD data define a rather 
steep spectrum with P = 3 — 4. 

The acquired HXD-PIN data are fully consistent with 
them coming from the same sources as detected with the 
XIS. Regardless of the way of removing the estimated con- 
tamination from X-2 and the fainter sources (section 3.4), 
the overall Suzaku spectrum, spanning 3-20 keV (with a 
narrow gap over 10-12 keV), is clearly more convex than 
a single PL, even though it can be approximated by a PL 
over a limited 3-8 keV range. This reconfirms previous 
broad-band studies with Ginga (Tsuru 1992), BeppoSAX 
(Cappi et al. 1999), and RXTE (Rephaeh, Gruber 2002), 
which all reported "thermal" type spectral shapes rather 



than straight PL like ones. Clearly, relativistic jet emis- 
sion, invoked by some authors (e.g., Krauss et al. 2005) 
to explain the behavior of some ULXs, is inappropriate at 
least in the case of M82 X-1, since the observed spectral 
curvature (|Ar| ^ 1 across ~ 10 keV) is too strong to be 
attributable to synchrotron cooling or burn-off effects. 

The curved spectrum of M82 X-1 can be expressed by 
either of the three successful continuum models; an em- 
pirical cutoff PL model, a thermal Comptonization model 
invoking a cool and thick electron cloud, and an optically- 
thick emission from a rather hot accretion disk (section 
3.4). Details of their model parameters depend to some 
extent on the assumed broad-band spectra of the contam- 
inants. Superposed on the continuum, we also detected 
two iron emission lines. 

4.2. An interpretation of the curved spectrum 

Although ULXs show a variety of spectral shapes, 
two prototypical spectral states have been recognized 
through past observations (e.g., Makishima et al. 2000; 
Sugiho 2003; Feng, Kaaret 2005; Dewangan et al. 2006; 
Stobbart et al. 2006; Mizuno et al. 2007; Makishima 
2007a; Makishima 2007b; Miyawaki 2008). One is "PL 
state" , wherein the spectrum below ^ 10 keV is approx- 
imated by a single PL model, but is often accompanied 
by a soft excess (e.g.. Miller et al. 2003; Soria et al. 2007) 
as well as a mild turn-over toward higher energies (e.g., 
> 6 keV; Mizuno et al. 2007). We already argued in 
section 3.4.3 that the second brightest source, M82 X- 
2, is likely to be an object in the PL state. The other 
is "convex-spectrum (hereafter CS) state" (Makishima 
2007a), in which the spectrum exhibits a significantly 
curved shape even in the limited energy range of ~ 1 to 

10 keV, so that an MCD model gives a better fit than a 
PL. The two ULX states are not considered to represent 
distinct classes of objects, since an increasing number of 
sources are confirmed to make transitions between these 
two spectral states (e.g., Makishima et al. 2000; Kubota 
et al. 2001b; Sugiho 2003; La Parola et al. 2004; Miyawaki 
2008). 

Mizuno et al. (2007) and Makishima (2007a) utilized 
Suzaku data on NGC 1313, to study spectral properties 
of two ULXs therein, X-1 and X-2. Then, X-2 was found in 
the CS state, and its 0.5-10 keV spectrum was represented 
successfully by a variable-p model with p = 0.63 ±0.03 and 
Tin = 1.86 ±0.15 keV. In contrast, the other source, X-1, 
was in the PL state, because its spectrum was much less 
convex (approximated by a PL with P ~ 2.1), with a hint 
of soft excess and a mild high-energy turn over. Thus, the 
XIS spectrum of NGC 1313 X-1 was reproduced success- 
fully by a cutoff PL continuum with Scut = 3.4lto4o keV 
and r = 0.89 ±0.20, plus a cool (Tin ~ 0.2 keV) diskbb for 
the soft excess. The source exhibited a rather high 0.4-10 
keV luminosity of 3.3 x 10''° erg s~^. 

As clarified in the present study, the Suzaku spectrum 
of M82 X-1 is approximated by a PL model with F ~ 1.7 
in the limited 3-8 keV range, while it is reproduced by 
a cutoff PL model with £;cut ^ 4 keV and F = 0.2 - 0.3 
(table 3) when the broad 3-20 keV band is incorporated. 
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These spectral properties of M82 X-1, as well as its lu- 
minosity, are very similar to those of NGC 1313 X-1 as 
recorded with Suzaku, although the values of T differ by 
~ 0.5. Based on these similarities, we conclude that M82 
X-1 was in the PL state during the present observations. 
This also implies that M82 X-1 can be regarded as an "or- 
dinary" ULX, rather than a more enigmatic class of ac- 
creting source, even though it is rather extreme in terms 
of the luminosity. 

While the cutoff PL is a purely empirical model with- 
out physical meaning, it approximates the process of un- 
saturated thermal Comptonization (section 3.3.3). In 
fact, the 3-20 kcV spectrum of M82 X-1 has also been 
reproduced reasonably well by the compTT model (ta- 
ble 3), invoking a rather low electron temperature of 
~ 2.3 keV and a large optical depth of ^ 10; the con- 
dition is close to a nearly saturated Comptonization. As 
mentioned in section 3.3.3, the parameters derived here 
agree very well with those of the same object obtained 
by Agrawal, Misra (2006) using the 3-10 keV XMM- 
Newton EPIC data. The Comptonization modeling has 
also provided reasonable account for the spectra, typi- 
cally in the 1-10 keV bandbass, of other PL-state ULXs 
(Kubota ct al. 2001a; Fiorito, Titarchuk 2004; Makishima 
2007a; Makishima 2007b; Soria ct al. 2007; Miyawaki 
2008). 

A valuable guideline to the understanding of ULXs may 
be provided by the study of Galactic black-hole binaries 
(BHBs) (Makishima 2007a). So far, a fair number of 
BHBs have been observed to leave the canonical standard- 
disk state, and enter "Very High" state (Miyamoto et al. 
1991; McClintock, Remillard 2006), when their luminosi- 
ties approach the corresponding Eddington limits. Then, 
the spectrum becomes dominated by a PL-like continuum 
with r > 2.4, plus a soft excess attributable to the un- 
derlying disk emission. Such examples include GX339— 4 
(Miyamoto et al. 1991), GS 1124-68 (Miyamoto et al. 
1993), GRO J1655-40 (Kubota et al. 2001a; Kobayashi 
et al. 2003), XTE J1550-564 (Kubota, Done 2004), 
GRS 1915 105 (Reig et al. 2003; Done et al. 2004), and 
4U 1630-47 (Abe ct al. 2005). These Very-High state 
spectra of BHBs have been interpreted successfully by 
assuming that a cool standard accretion disk becomes 
surrounded by a cloud of hot electrons, or a "corona", 
so that a significant fraction of the disk photons get 
Comptonized by the corona (e.g.. Done et al. 2007 and 
references therein). According to Kubota et al. (2001a), 
the RXTE PGA spectrum of GRO J1655-40 in the Very 
High state can be understood by invoking a corona with 
an electron temperature of To ~ 10 keV and a Compton op- 
tical depth T ~ 2, located above a cool disk which supplies 
the seed photons. Kobayashi et al. (2003) incorporated 
the RXTE HEXTE data of the same object and derived 
To = 75 - 85 keV and r = 0.45 - 0.65. Although the solu- 
tions thus degenerate, the Compton j/-parameter becomes 
in either case y ^ 0.3. 

After the initial suggestion by Kubota, Done, & 
Makishima (2002) and subsequent discussion by sev- 
eral authors (Tsunoda et al. 2006; Mizuno et al. 2007; 



Makishima 2007a; Miyawaki 2008), we presume that the 
PL state of ULXs is analogous to the Very High state 
of BHBs, and that both appear at near the Eddington 
luminosity. The bolometric luminosity of M82 X-1 dur- 
ing the present observation is estimated to be Lboi = 
(1.5 — 3) X lO'"' erg s~^, by integrating the cutoffPL model. 
The result does not differ significantly if we instead use 
the MCD fit. Then, its identification with the Eddington 
limit leads to a mass estimate as 100 — 200 Mq. The 
mass could be even higher by a few times, if considering 
that the highest recorded luminosity reaches ~ 1 x 10*^ 
erg s^^, and that the Very High state may appear at a lu- 
minosity somewhat below the Eddington value (Kubota, 
Makishima 2004). A possible formation scenario of such 
an intermediate-mass BH, particularly in close vicinities 
of galaxy nuclei, has been proposed by Ebisuzaki et al. 
(2001) and Portegies Zwart et al. (2006). We admit, on 
the other hand, that the reality of such a cool and thick 
Compton cloud needs a careful consideration, from the 
viewpoints of its stability, energetics, and geometry with 
respect to the seed photon source. These issues are dis- 
cussed by Agrawal, Misra (2006). 

4-3. Another interpretation of the spectrum 

Assuming that the PL state of ULXs is analogous to 
the Very High state of BHBs, we may identify, after 
some previous works (Mizuno et al. 2001; Watarai et al. 
2001; Tsunoda et al. 2006; Makishima 2007a), the CS 
state of ULXs to the slim-disk state of BHBs that is ob- 
served at their highest luminosities (Kubota, Makishima 
2004; Abe et al. 2005; Makishima 2007b). This analogy 
has the following three grounds. First, a ULX generally 
becomes more luminous, by a factor of 1.5-3 (in the 1- 
10 kcV band), when it makes a transition from the PL 
state to the CS state (Kubota et al. 2001b; La Parola et 
al. 2004; Sugiho 2003; Miyawaki 2008). This luminosity 
relation, including its sense and magnitude, is similar to 
that observed between the Very High state and the slim- 
disk state of BHBs. Second, the CS-state spectra of ULXs 
can be generally reproduced successfully by the variable-p 
model with p ~ 0.6 (e.g., Tsunoda et al. 2006; Mizuno et 
al. 2007). Finally, when multiple observations of a sin- 
gle CS-state ULX are analyzed by a simple MCD model, 
we observe the luminosity to scale as oc (Mizuno et al. 
2001). The second and third properties agree with the ob- 
served behavior of BHBs in their slim-disk state (Kubota, 
Makishima 2004), as well as predictions by the slim-disk 
theory (e.g., Watarai ct al. 2001; Kawaguchi 2003). 

In practice, the above state classification scheme of 
ULXs can be sometimes rather ambiguous, because a 
variable-p spectrum with a very high disk temperature 
takes a rather similar shape to a cutoff PL model. For 
example, Vierdayanti et al. (2006) applied the variable-p 
modeling to several PL-state ULXs observed with XMM- 
Newton. The present M82 X-1 data are subject to the 
same ambiguity, because the cutoff PL model and the 
variable-p model have been similarly successful. In fact, 
Okajima et al. (2006) interpreted an XMM-Newton spec- 
trum of this object (plus some other point sources in- 
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elusive) as arising from an extreme slim disk with a 
highly super-Eddington luminosity, formed around a BH 
of 20 — 30 Mq. Even utilizing the broader energy band 
realized by Suzaku, the degeneracy was not solved, and it 
remained after removing the source contamination. 

Let us briefly examine the slim-disk interpretation of 
the present Suzaku data for its consistency. The variable- 
p fit to the XIS plus HXD data gave extremely high disk 
temperatures as Tin = 3.0 — 3.5 keV (table 2 and table 3), 
regardless of the treatment of the source contamination. 
It significantly exceeds those observed from other PL-state 
ULXs, including 1.3-1.9 keV from NGC 1313 X-2 at a 
0.4-10 keV luminosity of 0.6 x 10^° erg s~^ (Mizuno et al. 
2007), and 1.5 - 1.7 keV from M81 X-9 at a 0.5-10 keV 
luminosity of 2.4 x 10^^° erg s^^ (Tsunoda et al. 2006). On 
the other hand, the temperature index of M82 X-1 was ob- 
tained as p = 0.64 ± 0.02 in table 2, or p = 0.75 in table 3. 
Since the latter is consistent with the condition of a stan- 
dard disk, we may calculate its innermost disk radius i?in 
via Lboi = 47r(i?in/0^f ^it; here a is the Stefan-Boltzmann 
constant, ^ ~ 1.2 is a correction factor in Makishima et al. 
(2000), and the disk inclination is assumed for simplicity 
as ~60°. 

Taking Ti„ = 3.0 keV (table-3) and Lboi = (1.5-3) x lO''" 
erg s~^ (section 4.2), the above procedure yields i?in = 
37 — 52 km, which corresponds to the last stable orbit of 
a non- rotating BH with a mass of only 4 — 6 Mq. Then, 
X-1 would be inferred to be shining at 25-35 times the 
Eddington luminosity. According to a recent numerical 
study by Vierdayanti et al. (2008), this mass estimate (and 
hence that of the super-Eddington factor) would not differ 
by more than a factor of 2, even if the system is in a slim- 
disk condition. Therefore, we must conclude the source to 
be in a very deep slim-disk condition. However, this con- 
clusion is not easily reconciled with the measured values 
of p, which is closer to the standard-disk condition of 0.75 
rather than to the limiting slim-disk condition of 0.5. As 
far as non-rotating BHs are considered, we hence regard 
the slim-disk interpretation as not being self-consistent; 
the observed spectrum is simply too hard to be regarded 
as emission from an optically-thick accretion disk. 

If the BH in M82 X-1 is spinning rapidly, its observed 
properties could be understood more easily in terms of 
standard accretion disks. This is because the last stable 
orbit would be then reduced down to ~ 1/6 times that 
of a non-spinning BH of the same mass. This idea was 
first proposed by Makishima et al. (2000) to explain CS- 
state spectra of ULXs, and recently applied successfully 
by Isobe et al. (2008) to a transient ULX. If M82 indeed 
hosts a maximally rotating BH, the BH mass estimated 
from the observed would increase to 24 — 36 Af©. 
However, we are still left with a super-Eddington factor 
of 4-9, and have to take into account various relativis- 
tic effects (Ebisawa et al. 2003). Therefore, the phys- 
ical reality of this interpretation is open. In addition, 
the reported spectral hardening toward lower fluxes (of 
presumably X-1 itself), based on the RXTE monitoring 
(Rephaeli, Gruber 2002), is opposite to the behavior of 
CS-state ULXs (Mizuno et al. 2007). From these consid- 



erations, we regard our baseline interpretation presented 
in section 4.2 as more appropriate. 

Incidentally, the state classification scheme of BHBs em- 
ployed above is based empirically on the observed lumi- 
nosity, but a more physically reasonable scheme should be 
based on the accretion rate normalized to the Eddington 
limit. Considering that the radiative efficiency itself 
varies, it is not obvious which of the two states, the 
Very High state and the slim-disk state, correspond to 
higher normalized accretion rates. They might even have 
nearly the same accretion rate, and are controlled by 
some other parameters including hysteresis effects. In 
GRO J1655-40 (Kubota, Makishima 2004), in fact, the 
system was observed to evolve directly from the standard- 
disk state into the slim-disk state on one occasion when 
the luminosity was increasing, while from the slim-disk 
state to the Very High state and then back to the 
standard-disk state on two occasions during the luminos- 
ity descent. A similar hysteresis effect, though less clear, 
is suggested in the case of 4U 1630-47 (Abe et al. 2005). 
Clearly, we need deeper understanding of these two high- 
luminosity states that have been recognized relatively re- 
cently. 

4-. 4- Subsidiary Spectral Features 

4-.4-1- Possible iron edge 

Even with the successful three continuum models found 
in section 3.3.3, we observe some data deficits in the 7 — 9 
keV energy range (figure 8). To account for this structure, 
we multiplied an absorption edge factor to the continuum 
models fitted to the XIS and HXD spectra. Regardless of 
the choice among the three successful continua, the value 
of then decreases by ~ 15 for a decrease of degree of 
freedom by 2, implying that the edge factor improves the 
fit with > 99% significance. The edge energy was con- 
strained as 7.1 ±0.1 keV, which agrees with the K-edge 
energy of nearly neutral iron. Its optical depth was ob- 
tained as 0.06 — 0.13, depending on the continuum models 
and the NXB systematic errors. 

The suggested excess Fe-K edge could be reconciled 
with the adopted low-energy absorption (1.1 x 10^^ cm~^) 
if the iron to alpha-clement abundance ratio were 5 to 10 
times higher than the solar ratio. However, this is oppo- 
site to the measured abundance pattern of X-ray emitting 
materials in M82 (Ranalli et al. 2008), and hence unlikely. 
An alternative possibility is that X-1 is in fact rather heav- 
ily absorbed by a column up to (4.5 — 9) x 10^^ cm~^ (as- 
suming solar ratios), but the associated strong low-energy 
absorption is masked by separate softer components. 

Yet another possibility is that the spectrum bears a 
"reflection" component from cold materials, and the Fe-K 
edge feature imprinted in that component is being ob- 
served. To examine this scenario, we removed the Fe-K 
edge factor from the cutoff PL continuum, and added to 
it a reflection component, pexrav (Magdziarz, Zdziarski 
1995), with the reflector iron abundance fixed at 0.6 so- 
lar, those of lighter elements at 1.2 solar (Ranalli et al. 
2008), and the inclination at 60°. We used the cutoff 
PL continuum itself as the input to the reflection. Then, 
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the fit chi-squared decreased (from the original cutoff PL 
modeling) by ^ 14, for a decrease of degree of freedom 
by 1. Therefore, the improvement is comparable to that 
obtained by applying the Fc-K edge. The reflector solid 
angle O, as seen from the nucleus, is suggested to be rel- 
atively large as ri/27r ~ 1, but with a considerable un- 
certainty, as the associated Compton hump is not well 
constrained by the HXD data. We consider this scenario 
most likely to be the case. 
4-4-^- Iron emission lines 

Let us also discuss the two Fe-K lines detected with 
the XIS (figure 7). Using Chandra and XMM-Newton, 
Strickland, Heckman (2007) showed that both these lines 
are associated primarily to diffuse emission, with insignif- 
icant contributions from bright point sources including X- 
1. Over a circular region of radius 29" (0.53 kpc) around 
the nucleus, they measured the 6.7 keV and 6.4 keV com- 
ponents to have photon fluxes of (6.9 ± 4.6) x 10~^ ph 
cm~^ s~^ and (4.8±4.0) x 10~^ ph cm~^ s~^, respectively 
(from their table 3). 

While our data integration radius (3') is significantly 
larger than that of Strickland, Heckman (2007), diffuse 
emission in the 2-8 keV band is mostly contained within 
the 30" radius (figure 4f of Strickland, Heckman 2007). 
Therefore, the Chandra and Suzaku results may be di- 
rectly compared to the first approximation. From our 
analysis conducted in section 3.3.1, the stronger 6.7 keV 
line has a fiux of 10.2 ±2.1 while the weaker 6.4 keV line 
4.5 ±2.0, both in units of 10"^ ph cm^^ s-\ Thus, the 
Suzaku and Chandra results agree with each other within 
respective errors. The observed Fe-K line photons are 
hence attributed primarily to the diffuse component, but 
some fraction of the 6.4 kcV photons may be associated 
with the putative refiection component of X-1. 

Intriguingly, the co-existing 6.4 kcV and 6.7 kcV Fc-K 
lines have also been detected in the diffuse X-ray emis- 
sion around our Galactic center (Koyama et al. 2007, and 
references therein) and along the Galactic ridge (Ebisawa 
et al. 2008). Therefore, we may be observing a similar 
phenomenon from the M82 central region. However, the 
emission in our Galactic center is confined within a con- 
siderably narrower region of ^ 150 pc, and its Fe-K photon 
luminosity, as measured with Suzaku, is 3 orders of mag- 
nitude lower than that from the central region of M82 
(Tsuru, Nobukawa, and Koyama, a private communica- 
tion). Presumably, the diffuse emission in M82 is signifi- 
cantly enhanced due to the extreme star-formation activ- 
ity. Further discussion on the iron emission lines is beyond 
the scope of the present paper. 

4-.5. Interpretations Invoking Anisotropic Emission 

In addition to the interpretations of ULXs assuming ei- 
ther intermeddate-mass BHs (section 4.2) or highly super- 
Eddington luminosities (section 4.3), a third explanation 
of these sources invokes strong emission anisotropy en- 
hanced toward us (e.g.. King 2002). Among such ideas, we 
have already argued in section 4.1 against the relativistic 
beaming scenario. The remaining possibility is geometri- 
cal coUimation of radiation by a funnel of thick materials 



(King 2002). In such cases, we should expect strong spec- 
tral signatures due to radiation reprocessing, including Fe- 
K edge features and fluorescent Fe-K lines, like in some 
extreme Seyfert 2 galaxies (e.g., Ueda et al. 2007). Indeed, 
we found the evidence of reflection by nearly neutral ma- 
terials (section 4.4). However, the observed reprocessor 
solid angle, fJ/27r'^ 1, and the equivalent width of the neu- 
tral Fe-K emission line, ~ 20 eV (section 3.3.1), are both 
far smaller than would be expected {n/2Tr >> 1, and an 
equivalent width of > 1 kcV) when the observed radiation 
is dominated by the reflected component. Furthermore, 
no adequate explanation based on the coUimation scenario 
has been presented, either for the two types of ULX spec- 
tra, or for the transitions between them. Therefore, the 
anisotropic radiation interpretation is unlikely to give a 
successful account of the ULX phenomenon in general. 

5. Conclusion 

Based on our discussion, we conclude that M82 X-1 was 
in the PL state during the present Suzaku observation, 
and was radiating roughly close to, or slightly below, the 
Eddington luminosity. This leads to a mass estimate of at 
least 100 — 200 Mq, which is in a good agreement with 
an independent estimate based on the QPO frequency 
(section 1). We therefore consider M82 X-1 as a genuine 
intermediate-mass BH, in the sense that it is significantly 
more massive than those BHs that can be explained by 
the standard stellar evolutionary scenario. 
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Table 1. Signal and background count rates from HXD-PIN. 





October 4 


October 19 


October 27 Average 


On-source 


0.295 ±0.006 


0.325 ±0.005 


0.343 ± 0.006 — 


NXB* 


0.260 


0.280 


0.300 — 


CXBt 


0.018 


0.018 


0.018 — 


Net signal^ 


0.016 ±0.006 


0.027 ±0.005 


0.025 ±0.006 0.023 ±0.003 



* Modeled with bgd.A. 

■f Predicted from the HEAO-1 measurements (Boldt 1987; Gruber et al. 1999). 
Net count rate, obtained by subtracting the NXB and CXB from the on-source rate. 



Table 2. Best-fit model spectral parameters with 90% confidence errors.* 



Models 


Parameters 


XIS-onlyT 


XIS+PIN^ 


NXB+4%!! 


NXB-4%i' 


Power-law 


r 


1.67 ±0.02 


1.75 ±0.02 


1.77±0.02 


1.73 ±0.02 




xVd.o.f 


318.6/174 


606.8/182 


809.1/182 


464.9/182 


Bremss 


T (keV) 
xVd.o.f 


14.81^7 
223.9/174 


292.2/182 


13.6t[!« 
408.3/182 


13.51^5 
242.0/182 


Cutoffpl 


r 

E,ut (keV) 
xVd.o.f 


L*- 'O-0.15 
5 8+1-0 

183.7/173 


„ __+0.08 
' -0.15 
c 7+0.5 
3- '-0.7 

193.6/181 


0.55in? 
228.6/181 


202.1/181 


CompTT 


Te (keV) 


2.5±0.1 


2.5±0.1 


2.4±0.1 


2.8±0.1 




T 


8.3 ±0.4 


8.0 ±0.4 


8.6 ±0.4 


7.4 ±0.4 




xVd.o.f 


184.5/173 


197.2/181 


203.5/181 


239.4/181 


Diskbb 


Tin (keV) 
xVd.o.f 


2.77 ±0.04 


2.79 ±0.04 


2.75 ±0.04 


2.80 ±0.04 




224.8/174 


243.4/182 


239.8/182 


312.8/182 


Variable-p 


Tin (keV) 
P 


q qQ+U.vi2 
0.64 ±0.02 


o c:o+U.28 

0.63 ±0.02 


q 7+0.22 
"J-*' -0.07 

0.63±0.02 


q cq+0.50 

0.63 ±0.02 




xVd.o.f 


183.4/173 


191.9/181 


225.1/181 


220.0/181 



* The absorbing column density is fixed at 1.1 x 10^^ cm ^. 
^ Using the XIS FI and BI spectra in the 3.2-10 keV band. 
i From the XIS and HXD-PIN joint fit. 

§ When the modeled NXB spectra are systematically over- or under-estimated by 4%. 
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Table 3. Same as table 2, but considering contamination from surrounding 



Models 


Parameters 


Case 1* 


Case 2T 


Power-law 


r 


1.64 ±0.03 


1.71 ±0.04 




xVd.o.f 


538.0/182 


303.9/182 


Bremss 


T (keV) 


18.2±1.1 


13.5 Q g 




xVd.o.f 


384.4/182 


239.0/182 


Cutoffpl 


r 

Ecnt (keV) 


0.27+lj-U^ 

— U. U4 

4.2 ±0.1 


0.17+^4 

— U.U4 

3.7±0.1 




xVd.o.f 


193.2/181 


193.1/181 


CompTT 


Tc (kcV) 
r 

xVd.o.f 


2.30 ±0.03 
9.7±0.2 
196.8/181 


9 qq+U.U4 

Z.OO Q Qg 

9 2+° '2 

—0 3 

199.9/181 


Diskbb 


Tin (keV) 
xVd.o.f 


3.09 ±0.07 
192.7/182 


^•■^3-0.10 

195.0/182 


Variablc-p 


Tin (keV) 
P 

xVd.o.f 


3.09 ±0.07 

7f^+0-0i 

^■'3-0.02 

192.7/181 


9 q^+O.UV 

^■^"3-0.11 
'3-0.03 

194.9/181 


Flux ^ 


(2-10 keV) 


1.2 


0.7 



* The case when spectra of the fainter contaminant sources (except X-2) were removed. 

The case when spectra of both the fainter sources and X-2 were removed. 
^ In units of 10"^^ erg s~^ cm~^. Depends no more than 10% on the spectral models. 
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Fig. 1. A background-inclusive 2-10 keV XIS image from the second observation. The smaller and larger squares show the 
FWHM and full-width at zero transmission FOVs of HXD-PIN, respectively. Dashed circles show the source and background regions 
used for the XIS analysis. A filled small circle indicates the location of the M81 nucleus, whereas two gray circles those of the two 
ULXs in this field, M81 X-6 and Holmberg IX X-1 (Fabbiano 1988a). 
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Fig. 2. XIS light curves in the 2-10 keV band extracted from the source (black) and background (gray) regions, acquired on (a) 
October 4, (b) October 19, and (c) October 27. The bin width is 256 s. 




Fig. 3. Background-subtracted XIS spectra of the core region of M82, averaged over the three observations. The instrumental 
response is included. The FI (three cameras summed) and BI CCD data are shown in black and gray, respectively. 
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Fig. 4. Earth-occultation spectra of HXD-PIN (black) obtained on (a) October 19 and (b) October 27, compared with the modeled 
non X-ray background (gray) to be observed during the same time periods. 
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Fig. 5. The 12-20 keV light curves of the M82 region acquired with HXD-PIN on (a) October 4, (b) October 19, and (c) October 
27. The observed raw data are shown in black, while the modeled NXB data in red, both using a binning of 256 s. Green data 
points show the signals remaining after subtracting the NXB and the expected CXB, binned into 1024 s. 
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Fig. 6. Time-averaged HXD-PIN spectra obtained on (a) October 4, (b) October 19, and (c) October 27. Each panel shows the 
observed spectrum (black), the modeled NXB spectrum (red), the signals remaining after subtracting the NXB (green), and the 
expected CXB spectrum (blue). 




Fig. 7. (a) A 3.2-10 kcV portion of the XIS FI (black) and HI (red) spectra in figure 3, presented together with the best-fit absorbed 
power-law plus gaussian model (upper panel) and the fit residuals (lower panel), (b) The same as panel a, but the power-law is 
replaced by a cutoff power-law. 
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Fig. 8. The same as figure 7, but the fit jointly incorporates the HXD-PIN spectrum (green) from which the NXB (nominal) and 
CXB are both subtracted. The employed fitting models are; (a) power-law, (b) bremsstrahlung, (c) cutoff power-law, (d) compTT, 
(e) diskbb, and (f) variable-p disk model. 
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Fig. 9. The same joint XIS plus HXD fits as figure 8, but incorporating estimated contributions from the contaminating sources 
(table 3). The data points from XIS-FI and HXD-PIN are shown in brown and green, respectively, while the XIS-BI data are omitted 
for clarity. The two Gaussians, though included, are omitted for the same reason, (a) The case when the summed spectrum of "the 
other fainter sources", estimated in the XIS band (blue), extends into the HXD-PIN range with a PL shape, (b) The case (Case 1) 
when the spectrum of the fainter sources (blue) is approximated by a 10 keV bremss model, with its normalization fixed in the XIS 
range. The X-1 spectrum (red) is modeled by a cutoff PL with a free shape and free normalization. The total model prediction is 
indicated in black, and the fit residuals are shown in the lower panel, (c) The same as panel b, but when the X-2 contribution is 
included as a fixed model (orange) represented by another cutoff PL model (Case 2; see text). In the HXD-PIN range, the three 
constituent components nearly overlap. 



